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Abstract. We prove that the smooth solutions to the Cauchy problem for the compressible 
general three-dimensional Ericksen-Leslie system modeling nematic liquid crystal flow with conserved 
mass, linear momentum, and dissipating total energy, generally lose classical smoothness within a 
finite time. 
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1. System, known results, and main problem. The liquid crystal state is 
often viewed as an intermediate state between liquid and solid. The molecules possess 
none or partial positional order but display a preferred orientation. The nematic phase 
is the simplest among all liquid crystal phases and is closest to the liquid state. The 
molecules float around as in a liquid, but have the tendency to align along a preferred 
direction due to their orientation. The hydrodynamic theory of liquid crystals, due to 
Ericksen and Leslie, was developed in the 1960's [H HI [211 [22] . The full Ericksen-Leslie 
system consists of the following equations (cf. [101 (Sll [23l [27l 119]): 

5tp -I- div(pu) = 0, (1.1) 
/ dW \ 

dtipu) +Div{pu(E)u) + Wp = -dj [p-Qj-''^d.j +Div{a), (1.2) 

Jp^^^ = (h- g) X d, u} = dx^d, (1.3) 

p = Ap'', yl>0, 7>1. (1.4) 

Equations (jl.l[) - (|1.4I) are given in M x R^^ and represent the conservation of mass, 
linear momentum, and angular momentum, respectively, with the anisotropic feature 
of liquid crystal materials exhibited in (|1.3I) and its nonlinear coupling in (|1.2p (cf. 

[21 [iU). 

Here, we consider the flow of a compressible isentropic material: p is the fluid 
density, p is a pressure, J is a positive inertial constant, 7 is an adiabatic exponent, 
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u — {ui, U2, Ws)"^ is the velocity, d = {di, d2, d^)"^ is the orientational order parameter 
representing the macroscopic average of the molecular directors, and := ^ + (w, V) 
denotes the material derivative. 

The notations A = i(VM + (Vw)"^), il = i(Vu-(Vu)^), N = ^d- fid represent 
the rate of strain tensor, the skew- symmetric part of the strain rate (the vorticity of 
the fluid), the material derivative of d (transport of center of mass), and the rigid 
rotation part of the changing rate of the director by the fluid vorticity, respectively. 

The Oseen-Zocher-Frank free energy functional W for the equilibrium configura- 
tion of a unit director field for a nematic crystal is given by the sum of the splay, the 
bend, and the twist, i.e, 

W ^ Ki^idiv df + K2^\d X {curl d)\^ + Ks^id ■ curl df > 0, (1.5) 



and the vector field 



is the molecular field. 

The kinematic transport g of the director d is defined by: 

g, = \iNi + X2djAji, (1.7) 

and represents the effect of the macroscopic flow field on the microscopic structure. 
The material coefficients Ai and A2 reflect the molecular shape (Jeffrey's orbit) and 
how slippery the particles are in the fluid. 

The stress tensor a has the following form: 

CTij = fiidkAkpdpdidj + fi2Nidj + fisd^N^ + ^iA^-j + ^ir^A^kdudj -f fied^Ajkdk- (1.8) 

The independent coefficients fix, ...,/i6i which may depend on the material and tem- 
perature, are called Leslie coefficients. 

The following relations are frequently introduced in the literature (cf. [22l I23j ): 



Ai2 + Ai3 = A^e - M5- (1-10) 

The identities (|1.9p are necessary conditions to satisfy the equation of motion iden- 
tically (cf. [ig. Section 6). The identity pTTUI is called Parodi's relation (cf. [5^), 
which is derived from Onsager's reciprocal relations expressing the equality of certain 
relations between flows and forces in thermodynamic systems out of equilibrium (cf. 
[3T]). Under the assumption of Parodi's relation, we see that the dynamics of an 
incompressible nematic liquid crystal flow involve flve independent Leslie coefficients. 

As was proved in [TO], \d{0)\ = 1 imphes \d{t)\ = 1 and (c?(0),w(0)) = implies 
{d{t),Lj{t)) = for all time. 

When c? is a constant vector field, the system (|l.ip - (|1.4p becomes the compress- 
ible Navier-Stokes equations describing the motion of a compressible fluid. 

Since the mathematical structure of the Ericksen-Leslie system is quite compli- 
cated, the study of the full model presents several mathematical difficulties. Almost 
all existing investigations were restricted to simplified versions. For the incompressible 
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model in [24 , Lin introduced a simplification of the general Ericksen-Leslie system 
that keeps many of the mathematical difficulties of the original system by using a 
Ginzburg-Landau approximation to relax the nonlinear constraint d G 5^. Namely, 
instead of the restriction \d\ = 1 he added the penalty term — 1)^ to the free 

energy functional W and neglected the inertial constant J. Later in |25J , Lin and Liu 
showed the global existence of weak solutions and smooth solutions for that approx- 
imation. In [35j , a very simple proof of local well-posedness for this coupled system 
was provided using a contraction mapping argument. It was proved that this system 
is globally well-posed and has compact global attractors in 2D. For more results on 
the approximative system see [26l |28l |30l |36l [HI H]- Further, taking the limit of 
e in the penalty term, one can get a coupling between the compressible Navier- 
Stokes equations and a transported heat flow of harmonic maps into [27j . It is 
a macroscopic continuum description of the evolution for liquid crystals of nematic 
type under the influence of both the flow field u, and the macroscopic description 
of the microscopic orientation configurations d of (rod-like) liquid crystals. Recently, 
Hong 14 and Lin-Liu- Wang [29 showed independently the global existence of weak 
solution of an incompressible model in two dimensional space. Moreover, in [29], the 
regularity of solutions, except for a countable set of singularities whose projection on 
the time axis is a finite set, has been obtained (see also [H]). In [37], Wang established 
a global well-posedness theory for the incompressible liquid crystals for rough initial 
data, provided that HwoUbmo-i + II'^o||sa/o < eo for some eo > 0. In [4], regularity 
and uniqueness for solutions to density dependent nematic liquid crystals systems in 
the Ginzburg-Landau approximation were established for a bounded domain: in 2D 
the system has a global classical solution, in contrast to the 3D case. A family of exact 
solutions with finite energy to the incompressible liquid crystals in two dimensions 
was constructed in [7]. 

Concerning the compressible 3D case, local existence and uniqueness of strong 
solutions for the coupling between the compressible Navier-Stokes equations and a 
transported heat fiow of harmonic maps was proved (see jl8p. provided that the initial 
data po,uo, do are sufficiently regular and satisfy a natural compatibility condition. A 
criterion for possible breakdown of such a local strong solution at finite time was given 
in terms of blow up of the L°°-norms of p and Vd. Alternative blow-up criteria were 
derived in terms of the L°°-norms of Vu and Vd in [IT] and in terms of integral of L°°- 
norms of Vu and the BMO-norm of Vd, in [3] . The global existence of weak solutions 
with large initial data is still an outstanding open problem for dimensions larger then 
or equal to 3. So far, only results in one space dimension have been obtained, for 
instance, we refer to [5J |S]. In [IB], the existence and uniqueness of global strong 
solutions to the Cauchy problem is proved for 3D in critical Besov spaces provided 
that the initial data is close to an equilibrium state. 

The overview of numerical methods used for the nematic liquid crystal flows can 
be found in [1 . 

At the same time, up to now there are no results on the classical solvability for 
the main initial and initial-boundary problems for the full system (|l.ip - (|1.4p . In the 
present paper we prove that a global classical solution to the Cauchy problem does 
not exists, and that this phenomenon is due to the presence of the viscosity term. The 
nature of loss of smoothness is similar to the case of the compressible Navier-Stokes 
equation and the anisotropic features do not influence this phenomenon. 
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P(t) = pudx 



2. A general blowup result. In this paper we study the system (|l.ip - (|1.4p 

with the initial condition 

{p,u,d) ={pa,uo,do), do G 5^. (2.1) 

t=o 

We introduce the following natural functionals defined on the solution of the 
system - (|1.4|) : mass 

■m{t) — / pdx, 



linear momentum 



and total energy 

E{t) = \ I p\u\^dx^ 5— ( pdx + \ ( Jp\uj\^dx+ [ pW dx = 

i;fc(t) + E,{t) + Ep^iit) + Ep^2{t) > 0, 

where Ek{t), Ei(t), and Ep^i, i = 1,2, are the kinetic, internal, and potential compo- 
nents of energy, respectively. 

Definition 2.1. A solution {p,u,d) to the Cauchy problem hl.l]] - |j.^[ ), V2. 1]) 
belongs to the class R if the solution is classical, p > 0, the mass m{t), linear mo- 
mentum P{t), and total energy are finite for all t > 0, and, in addition, the mass and 
linear momentum are conserved, i.e., m{t) — m — const, P{t) ~ P = const. 

Thus, if the solution belongs to the class ^, then 

p e n , ^fpu G L^, ^/puj e L^, Vrf e L^, {d, curld) G , d x curld G L\ 

where L'' = L'^(M^). 

Further, if we additionally assume 

ueH\ d^AdeL^, AdeL^, NgL"^. (2.3) 

A direct calculation with smooth solution {p, u, d) to the system (jl.ip - (|1.4p yields 
(cf. [17], Theorem 2.1), 

^^E{t)=- J (^p,\d^Ad\^ + ^\Vu\^ + {p, + pe)\Ad\^) dx 

R3 

+ Ai / |7Vp dx + (A2 -P2- P3) [ {N, Ad) dx. (2.4) 



R3 R3 

Here and throughout, we always assume that 

Ai < 0, P5+ M6 > 0, pi> 0, p4 > 0. (2.5) 
These assumptions are necessary conditions for the dissipation of the director field 

[ni[53]. 
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Indeed, as was shown in [S^ by means of the Holder inequahty, if (|2.5p holds, 
under the additional assumption 



IA2 - ifJ-2 + M3)| < 2y/-Xi{ll5 + He) 



(2.6) 



one has the following energy inequality: the basic energy law (without Parodi's rela- 
tion) holds 



d_ 



E{t)<-j (^i|d^Ad|2 + ^|VM|2) dx<0. 



(2.7) 



Our main result is the following theorem. 

Theorem 2.2. Suppose that 7 ^ | O'nd inequalities i2.5]) and \2. 6\) hold. Then 
there is no global in time solution to the Cauchy problem il.l]) - Ji.^p , i2.1\) in the 
class ^ satisfying 112. 3\) . 

The proof is similar to |33j , |34j , where a blow-up result was proved for the com- 
pressible Navier-Stokes system and the compressible magnetohydrodynamics system. 

The next lemma is a key technical point of the proof of the theorem. 

Lemma 2.3. Let 7 > | and u £ i?^(K^). // \P\ ^ 0, then there exists a positive 
constant C such that for the solutions from the class R the following inequality holds: 



iDul^dx > CE- 



(2.8) 



Proof. First, from the inequality 




u\— dx\ <Ci \Du\^dx 



(2.9) 



where the constant Ci > depends on n, n> 3, ([13J, p.22) we get H^W) C L^{R^). 
Then, using 7 > I, the Holder inequality gives 



IPI 



pu dx 



R3 




(2.10) 



Then, by the Jensen inequality (e.g., [20], Theorem 8.1.3), we have for 7 > |, 



5(7-1) 



J p'^ dx 



p 5 dx 



< 



mA 



(2.11) 



Thus, ([2l0| and ([CT]) imply 

|P| < C2 {E,{t))wrTj 




(2.12) 



with the constant C2 = mi (^) > 0. Thus, the inequality (|2?8l) follows from 

-cT 



T9\i . (|2T2l) . with the constant C = . □ 



Proof of Theorem [Ml Since Ei{t) < E{t) < E{0) by ([231), from (HH) and 
(|2.8p we conclude 



^E{t) < ^fiiC{E,{t))-wrTJ < -fiiC{E{0))-^. 
at 

which contradicts the non-negativity of E{t) for all t > 0. This concludes the proof 
of the theorem. □ 
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